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ABSTRACT 

Numerical  weather  analysis  relies  heavily  on  numerically  computed 
winds  both  for  direct  output  and  internally  for  calculating  advection  and 
other  purposes.     It  is  the  intent  of  this  paper  to  compute,  analyze,  and 
compare  the  geostrophic,  the  non-divergent,  the  divergent ,  and  a  com- 
bination of  the  latter  two.    Where  there  occurs  a  significant  difference 
between  these  numerically  derived  winds,  an  attempt  is  made  to  com- 
pare each  to  the  observed  winds.    This  study  is  made  for  all  standard 
levels  through  300  mb.     Data  and  several  computer  programs  were  fur- 
nished by  the  U»  S.  Navy  Fleet  Numerical  Weather  Facility,  Monterey, 
California . 

The  writer  wishes  to  express  his  appreciation     to  Professor  George 
Jo  Haltiner  of  the  U.  S.  Naval  Postgraduate  School  for  his  guidance, 
contributions  and  encouragement  in  this  work. 

Appreciation  is  expressed  to  Lieutenant  Commander  Mildred  J.  Seljos 
United  States  Navy  for  aid  and  advice  in  programming.    Further  thanks 
are  extended  to  Lt.  R.  F.  Alden,  U.  S,  Navy,  and  Lt„  G.  C.  Rosenber- 
ger  who  provided  programs  for  computation  of  horizontal  velocity  diver- 
gence and  for  advection  of  temperature. 


ii 


TABLE  OF  CONTENTS 

Section  Title  Page 

1 .  Introduction.  ............. 1 

2 .  Background  ....................  2 

3 .  Procedures 5 

4 .  Results 10 

5 .  Conclusions. 15 

6.  Illustrations 16 

7.  Bibliography „...o..»...  31 

8.  Appendix  I 32 


in 


Table  of  Symbols 
g  the  upward  component  of  the  apparent  gravitational  acceleration 

D  the  altimeter  correction       d  =  Z  -Z 


P 


Zp  the  height  of  an  isobaric  surface 


Z  the  pressure  altitude  in  the  standard  atomosphere 

f  the  Coriolis  parameter,  2   CI  sin0,  where     0    is  the  geographi- 

cal latitude 

f  the  mean  value  of  the  Coriolis  parameter,      2  Q,  sin  45° 

$  the  geopotential 

r  the  relative  vorticity         £=  A_v_  „  .liL 
s  dx  dy 

Tj  the  absolute  vorticity        <n  =  £+  f 

^  the  stream  function  for  the  non-divergent  component  of  velocity 

^  the  velocity  potential  for  the  divergent  component  of  velocity 

to  the  vertical  velocity  of  air  in  pressure  coordinates 

P  the  atmospheric  pressure  in  millibars 

,_.i.  T       5  6 

a  the  static  stability  a=  ~~      ^n 

R  the  universal  gas  constant 

d  grid  distance,  381  km  at  60N 

m  map  scale  factor 

v  del  operator  on  a  constant-pressure  surface 

y2  horizontal  Laplacian  operator  on  a  constant- pressure  surface 

J  horizontal  Jacobian  operator    J(ft;BJ"<-)x     <ty      ^y  «** 
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^  a  finite  difference  Laplacian  V^fain  ^/C+^i-i  +aT'-  4Ai  ) 

ir  finite  differnece  Jacobian  r>  J     aj  wJ+i  v'-M   /jt<  J-<)/Dj     dj  \ 

Vg  geostrophic  wind 

Ve  non-divergent  wind 

V^  Ve 

Vx  divergent  wind 

Vt  total  wind 

Va  actual  wind  as  measured  from  soundings 

£  dummy  variable 

I£  wind  speed  of  V£ 

DGD£  difference  in  wind  direction,  Dg  -  D£ 

IGl£  difference  in  wind  speed,  Ig  -  # 

DADV  advection  of  absolute  geostrophic  vorticity  by  Vg 

EADV  advection  of  absolute  stream  function  vorticity  by  Ve 

TADV  advection  of  absolute  stream  function  vorticity  by  Vt 

XADVT  advection  of  temperature  with  Vx 

t  time 

u£  u  component  of  V£ 

v£  v  component  of  V£ 

D£  wind  direction  of  V£ 
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1 .     Introduction 

One  of  the  primary  products  of  any  numerical  weather  facility  is 
wind  data  for  operational  use  and  for  further  numerical  computations 
such  as  advection.     It  is  the  purpose  of  this  paper  to  analyze  and 
compare  several  types  of  numerically  derived  winds  „    The  types  con- 
sidered are  the  geostrophic,  Vg,  the  non-divergent,  Ve,  the  divergent, 
Vx,  and  finally  the  total  wind,  Vt,  which  is  a  combination  of  the  di 
vergent  and  non-divergent  parts  of  the  wind.    These  winds  are  cal- 
culated by  solution  of  appropriate  equations  using  data  from  the  U.S. 
Naval  Fleet  Numerical  Weather  Facility,  Monterey,  California  (FNWF). 
Computations  were  made  on  the  U„  S.  Naval  Postgraduate  School  CDC 
1604  digital  computer. 

Comparison  of  these  winds  is  made  at  all  standard  levels  through 
the  300   mb  level  by  calculation  of  speed  and  geographical  direction  on 
separate  maps  for  each  level  covering  the    Northern  Hemisphere.    Dif- 
ference fields  of  direction  and  speed  are  made  between  each  wind  and 
the  geostrophic  wind.    Vertical  cross  sections  were  computed  and 
studied  for  each  of  the  winds  and  for  the  difference  fields.    Further  com- 
parisons and  analyses  are  made  by  using  these  three  winds  to  advect 
absolute  vorticity.    Advection  of  temperature  was  also  computed  using 
Vx. 


2 .     Background 

The  geostrophic  wind  is  computed  from  the  well  known  formula 
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Vj   =  "~  _f    \  dy  L  dx  &  >  (!) 

utilizing  a  standard  program  of  FNWF.    This  program  generates  fields 
of  geostrophic  wind  directions  in  degrees  r  id  speed  in  knots.  All  com- 
putations were  done  on  the  FNWF  square  grid  centered  on  the  North 
Pole  and  consisting  of  3968  with  a  grid  in  erval  of  381  km  at  60N. 
In  this  program  the  sine  of  latitude,  as  required  in  computation  of  the 
coriolis  parameter,  f,  was  set  equal  sine  15  degrees.    The  latter  points 
were  not  utilized  in  the  evaluation  of  results. 

Geostrophic  wind,  Vg,  is  utilized  as  a  basis  of  comparison     for 
later  derived  non-divergent „  Ve,  and  total,  Vt,  winds. 

The  Helm  holt-    theorem  allows  expression  of  the  wind  as  the  sum 
of  the  divergent  and  non-divergent  components. 

v=£xv^  +  vx.  (2) 

Winds  of  this  form  will  be  h   i   ■"■ ' '.-.ti  referred  to  as  total  wind,  Vt. 

The  potential  function  X  for  the  divergent  part  of  the  wind  is  obtain- 
ed  through  relaxation  of  the  Pois son-type  continuity  equation 

V2X  +  ff~-0.  0) 


The  stream  function,  >fc,  for  the  non-divergent  portion  of  the  wind  is 
approximated  by  solution  of  the  balance  equation 

Here,  $  is  the  geopotential. 

The  vertical  component  of  vorticity,  C  ,  is  the  Laplacian  of  the 
stream  function 

f=^  (5, 

and  absolute  vorticity,    r\    ,  is  the  sum  of  fand  the  vertical  component 
of  the  earths  vorticity,  fr 

(6) 


rj=  "f  +•£  ■ 


Vertical  velocity,  to,  is  computed  by  the  diagnostic  equation 

V^)tP^f^--}[vLJ(olr)-J(o^r)-fj(T,Ti)]  (7) 

developed  by  Haltiner,  Clarke,  and  Lawniczak   [l]. 

For  convenience  the  stream  function  has  been  redefined  so  as  to  be 
in  units  of  height, 

y-fE.  (8, 

The  geostrophic  advection  of  absolute  geostrophic  vorticity  is  com- 
puted from  the  expression 


and  the  advection  of  r\e  by 

The  advection  of  r\e  by  Vx  is  computed  from  the  formula 

-^2e     .     ,r    <*9?e  (11) 
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(10)  and  (11)  are  combined  to  represent  advection  of  r\e  by  Vt  and 
compared  to  similar  advection  using  Vg  and  Ve 


3       Procedures 

The  sequence  of  calculations  programmed  started  with  solution  of 
the  geostrophic  wind  equation   in    finite  difference  form 

(12) 


\V=  -,-Z7(D/-Di£-(D.-*. -l>>-i)ij 


ltd 

Fields  of  u    and  v    relative  to  map  projection  were  generated  and  convert- 
ed into  fields  of  wind  direction  in  geographical  headings  at  each  grid 
point.     Likewise,  fields  of  iso'achs  were  produced  and  printed  in  con 
toured  map  form.    These  calculations  and  maps  were  made  for  all  stand- 
ard levels  to  300  mb  for  this  and  all  subsequent  programs. 

The  balance  equation,   (4),  as  programmed    by  M„  Reese    of  FNWF, 
was  then  solved  and  the  fields  of  E  values  were  stored  on  magnetic 
tape.    The  iterative  scheme  of  this  program  does  not  always  converge  at 
300  mb  or  has  an  excessive  running  time,  both  of  which  limited  the  data 
produced.    The  E  values  were  then  used  as  input  to  the  geostrophic 
wind  program  with  the  appropriate  substitution  of  f  for  f.    These    "balanced" 
winds  were  computed  and    stored  in  the  same  manner  as  Vg. 

Divergence  fields  are  taken  as  input  for  solution  of  (3),     Here,  as 
with  >&,  X  is  redefined  as 

*  •  j  *  <i3> 


so  that  the  velocity  potential  has  units  of  height  and  therefore  is  com- 
patable  with  E  and  D  f or  u    sin  (12)  with  minor  changes.    Expressing 
the  first  term  of  equation  (3)  in  finite  difference  form  gives 

V  a-  ~     J   m2-     dp  (14) 

which  may  be  solved  by  use  of  an  FNWF  Poisson  relaxation  program. 
The  derivitive  -*d  tc/  d  p    which  is  the  horizontal  divergence  of  the  wind 
is  evaluated  by  a  finite  difference  approximation  as  programmed  by 
Alden  and  Rosenberger    [2J 

T1-*    £  fields  are  then  used  to  compute  the  divergent  wind  component 
by  the  finite  difference  equation 


ftL-*to+(sr-xif'#I        'i5> 


Again  the  u  and  v  components  are  stored  and  maps  of  direction  and  of 
isotachs  are  produced. 

Now  the  components  of  Ve  and  Vx  are  added  and  a  field  of  these  com- 
bined components  produced  as  follows 

u  \[r+  ux  =  ut  (16. a) 

v^+vx  =  vt  (16. b) 

As  before  these  component  fields  are  stored  and  maps  of  each  level  print 

ed. 

Programs  to  aid  comparison  and  analysis  are  now  written.  Analysis 
of  directional  differences  is  made  by  using  the  direction  fields  of  Vg  as 
a  basis  and  subtracting  directions  of  Ve  and  Vt  from  it.    This  proved  to 


be  a  challenging  problem  for  several  reasons.    First,  directions  are 
by  nature  discontinuous  at  singular  points  on  isogons  maps     and  also 
in  general  at  360  degrees.     Second,  direction  tends  to  vary  widely 
in  the  vicinity  of  small  scale  features  that  are  not  concentric  between 
the  various  fields  compared. 

The  discontinuity  problem  was  solved  successfully  with  a  rather 
intricate  program.     However,  the  singular  points  and  non-concentric 
small  scale  features  still  caused  a  large  gradient  of  direction  and 
difficult  contouring  to  plot  or  read.    This  problem  was  alleviated  by  re- 
placing a  value  of  minus  zero  at  all  points  on  the  directional  difference 
maps  at  which  the  wind  speed  was  less  than  five  knots.    This  improved 
appearance  of  the  maps  and  facilitated  evaluation „     Moreover  the  min- 
us zero  values  still  indicated  those  areas  where  singular  points  had 
occurred.    The  5-knot  criterion  proved  to  be  of  no  detriment  in  subsequent 
use. 

Speed-difference  fields,  IGIE  and  IGIT,  are  produced  by  the  same  pro- 
gram as  directional  difference.    The  computation  of  speed  differences  in= 
eluded  a  statistical  analysis  and  grid  maps  of  output  included  the  pillow, 
P,  and  root-mean-square  difference,  RMSD.     These  terms  are  defined  in 
APPENDIX  I. 

Another  series  of  programs  were  written  to  calculate  the  advection  of 
absolute  vorticity,  V  ,  by  the  various  winds  previously  computed  and  to 
compare  the  results. 


First,  fields  of  D  values  are  used  as  "input  to  evaluate  the  geostro 
phic  vorticity  -v)  3  {23     VT7  * 


^-"SfV'O  +  f  f17, 

with  an  FNWF  program.    Then  the  resulting  fields  from  (17)  are  used  to 
obtain 


a^Ve^- 


which  is  plotted  on  standard  grid  maps.    This  program  is  then  internal 
ly  modified  to  solve 

and 

(20) 


^  f-^fJJfE.n.) 


Next  previously  stored  components  of  Vx  are  used  to  compute  the 
advection  by  the  divergent  wind  as  follows: 

This  result  is  added  to  (2  0)  to  produce  fields  of  TADV. 

Next  these  three  advection  schemes  are  compared  using  DADV  as  a 
basis  and  forming  difference  fields  similar  to  those  used  for  direct  com- 
parison of  the  winds . 

Typical  computer  operating  time  including  efficient  handling  of  tapes 
and  printing  of  output  on  an  IBM  1402  printer  is  shown  in  Table  I. 
.    2s  shown  are  for  one  reporting  period  of  one  day  at  all  levels. 
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Table  I 


Program 


Run  time  in  minutes 


Solution  of  Balance  Equation,   (4) 

Wind  speed  and  directions  from  various  inputs 

Solution  of  equation  (14)  for  v 

Solution  of  (16. a)  and  (16. b) 

Comparison  of  winds 

Advection  of  absolute  vorticity  with  all  winds 

Comparison  of  advections 

Preparation  of  grid  field  for  FNWF  plotter 


46 

18  feach) 
30 
12 
35 
35 
20 
4  (each) 


Total 


200  minutes . 


It  is  apparent  that  operating  time  is  a  significant  factor  in  restricting  the 
amount  of  data  collected. 


4.    Results 

Results  and    conclusions  are  based  on  analyses  of  data  selected 
from  the  period  February  10  to  20,   1964,     Only  the  United  States  is 
used  in  figures  as  it  represents  an  area  of  reliable  data.     Using  a 
small  area  for  figures  also  allows  display  of  finer  detail.    For  illustra- 
tion 0000 Z  16  Feb  1964  was  selected  as  it  contains  a  wide  range  of 
synoptically  interesting  features  over  the  United  States. 

Figure  1  depicts  the  pressure  analysis  (dashed  lines)  superimposed 
on  the  stream  function  analysis  at  500  mb.    The  relative  height  differ- 
ence at  the  semi-enclosed  low  center  south  of  the  Great  Lakes  is  typ- 
ical of  the  500-mb  level  throughout  the  period  studied.     Figure  2  shows 
stream  function  analysis  of  the  same  period.     Comparison  with  pressure 
height  analysis  i  .dicate?  that  the  patterns  are  similarly  related  at  other 
standard  levels. 

Wind  speeds  computed  from  5  00-mb  height  and  stream  function  analy- 
sis are  shown  in  figure' 4 „    Dashed  isotachs  enclose  centers  of  maximum 
and  minimum  speeds.     Speeds  for  other  standard  levels  are  shown  in 
figures  2  and  3  where  dashed  lines  are  isotachs „     Large  numerals  indicate 
values  of  minimum  and  maximum  speeds.    Analysis  of  speeds  of  these 
two  winds  shows  there  is  a  significant  difference  between  them.     Dashed 
lines  in  figure  5  (lower)  are  isolines  of  difference  in  wind  speeds,  gao- 
strophic  minus  non-divergent.     Similar  analysis  of  other  levels  are  shown 
in  figures  6  and  7 . 
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Figure  5  (upper)  shows  differences  in  direction  for  500  mb,  Dg   -De. 
Solid  lines  show  the  500-mb  pressure  height  analysis,  alternate  dashed- 
dot  lines  are  contours  of  negative  5-kt  speeds,  and  dashed  lines  are  5-kt 
contours  of  positive  difference ,     The  area  marked    >  30  is  almost  com- 
pletely within  an  area  of  wind  speeds  of  less  than  5  kt.      Difference  in 
directions  are  not  shown  for  other  levels  but  are  easily  located  by  refer- 
ence to  figures  2  and  3  ,     Pillow  and  RMSD  values  for  the  hemispheric 
maps  is  given  in  figures  5  through  7.     Note  that  P  indicates  that  on  the 
hemispheric  scale  Ve  speed  is  lower  at  upper  levels,    e.g.  by  a  factor 
of  approximately  10  per  cent  at  300  mb.    Except  in  regions  of  very  sharp 
curvature  or  near  singular  points,  DGDE,  figure  5  (upper),  is  generally 
less  than  ten  degrees  at  this  and  all  other  levels.    A  more  detailed 
analysis  was  made  of  gridpoints  and  actual  winds  from  185  soundings 
in  areas  where  DGDE  was  greater  than  10  and  both  Ve  and  Vg  exceed 
5  kt.    The  graph  in  figure  8  shows  the  rslationship  found.     It  must  be 
pointed  out  that  interpolation  within  grid  points  was    required  to  use 
actual  soundings  and  that  actual  soundings  are  transmitted  to  the  nearest 
10  degrees  only.     In  an  effort  to  at  least  partially  reduce  effects  of  fric- 
tion and  local  terrain  effects  the  surface  winds  were  not  included  in  the 
data  used  for  construction  of  this  graph.     "Va"  is  the  actual  wind  as  taken 
from  the  soundings  and/ or  streamline  analyses.    A  similar  detailed  analysi; 
of  speed  differences  was  made  using  210  soundings  from  areas  where 
IGIE  exceeded  10  knots.    Results  are  shown  in  figure  9. 
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Maps  of  X  for  1000  mb  (lower)  and  300  mb  (upper)  are  shown  in  figure 
10 -     Dashed  lines  are  5  meter  contours,  solid  lines  are  2-kt  intervals 
of  wind  speed,   Ix.    Figure  11  shows  Xand  Ix  for  500  mb  and  850  mb. 
These  maps  were  plotted  using  an  FNWF  program  designed  for  another 
purpose  so  that  the  contour  labels  have  no  intrinsic  meaning.  Contour 
interval  is  10  meters.     Note  that  the  pre-trough  region  in  low  levels 
is  associated  with  a  \  maximum  (point  of  Vx  inflow)  at  upper  levels  and 
a  X  minimum  (point  of  Vx  outflow)  as  would  be  reflected  in  divergence 
and  convergence  patterns  associated  with  such  height  patterns. 

Ix  can  be  seen  to  be  maximum  at  the  300-mb  and  1000-mb  levels. 
Note  that  these  maximums  are  of  the  same  magnitude,  8  to  10  knots, 
but  that  directions  are  opposite.    As  might  be  expected,  speeds  are  min- 
imum at  500  mb    as  a  reflection  of  the  minimum  values  of  divergence 
known  to  exist  there. 

At  300  mb  even  maximum  values  represent  a  very  small  percent  of 
the  wind  speed,  seldom  over  10  percent  and  generally  less  than  three 
percent  of  Ve.    With  the  lower  wind  speeds    found  at  lower  levels  these 
same  magnitudes  of  Ix  are  found  in  some  instances  to  equal  45  percent 
of  Ve.     Such  a  case  is  depicted  along  the  Gulf  Coast  in  figure  8  (lower). 
It  is  of  further  interest  that  a  deepening,  fast-moving  cold  front  was 
located  in  the  same  area  almost  directly  under  the  1000-mb  Ix  maximum 
center.    This  front  moved  rapidly  up  the  East  Coast  and  into  the  North 
Atlantic o    A  center  of  maximum  Ix  was  associated  with  this  front  through- 
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out  the  traversal  of  the  United  States  but  speed  of  movement  and  the 
short  period  of  accurate  data  prevented  detailed  study  of  this  associa- 
tion o 

Combining  Ve  and  Vx,  Vt  was  formed ■     Figure  12  (upper)  shows  the 
effects  of  this  combination  at  850  mb;  12  (lower)  depicts     1000  mb 
Solid  lines  are  isotachs  of  Vt,  dashed  of  Ve„    The  effect  of  this 
combination  at  middle  and  upper  levels  is  so  small  as  to  be  undetectable 
in  this  type  analysis  and  is  not  shown.     The  effects  at  1000  mb  are  quite 
evident.     Detailed  analysis  of  these  proved  that  here  there  is  a  change 
in  location  of  relative  differences  but  no  change  in  magnitude  of  the 
difference.    Terrain  and  friction  effects  preclude  any  meaningful  com- 
parison  of  the  small  deviations  from  Ve  to  actual  wind  data.    Further 
attempts  to  contrast  Vt  with  Vg  only  further  emphasizes  that  there  is  no 
significant  difference  between  Vt  and  Ve, 

Figures  12  and  13  depict  the  results  of  using  Vg  and  Ve  for  advection 
of  absolute  vorticity.     Solid  lines  are  lines  of  zero  advection  using  Vg, 
DADV;  dashed  lines  are  for  Ve,  EADV.    The  300-mb  level  is  omitted  but 
showed  no  marked  difference  between  the  two  advecting  winds  as  it  too 
showed  practically  congruent  patterns.    The  numerals  in  the  ratios  shown 

are  located  at  points  of  local  minimum  and  maximum  value,  the  upper  num- 

=7  -1 

eral  is  for  DADV.     Units  are  10       sec       .     These  values  of  change  of 

vorticity  are  typical  of  the  period  investigated  and  showed  DADV  to  be 

larger  at  both  extremes  by  a  factor  of  5  to  10  units .     Computing  the  pillow 
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and  RMSD  for  all  the  standard  levels  gave  a  pillow  O  in  every  case  and 
RMSD  of  2  to  7  units  with  largest  being  at  300  mb„     Use  of  Vt  made 
no  detectable  difference  in  overall  patterns .     Some  minor  increases 
in  value  in  the  vicinity  of  sharp  curvature  at  low  levels  was  noted  but 
were  of  the  order  of  only  10  to  20  percent,    Advection  of  vorticity  with 
.x  was  not  computed  as  its  significance  has  been  seen  to  be  small  by 
its  small  effect  when  added  to  Ve.     However,  advection  of  temperature 
using  Vx  was  computed  using  a  program  which  evaluates  the  expressions 

J(E,T)  (21) 

and 

W-^T     .  (22) 

One  hour  change  of  temperature  was  computed „    Figure  14  shows  1000-= 
mb  change  of  temperature  with  Ve  (solid  lines)  ,   Units  are  degrees 
colsius  10™    .    Detailed  analysis  of  the  grid  map  shows  that  advection 
with  Vx  over  the  entire  map  was  less  than  *  0o01  degrees  over  the  United 
States    except  in  the  area  along  the  Gulf  Coast  where  values  are  depict- 
ed with  a  dashed  contour  enclosing  an  area  of  -  0,2  degrees .     Here  the 
advection  with  Vx  actually  exceeded  EADV,    This  is  the  area  noted  earl- 
ier of  the  deepening  low  pressure  system  and  front,,    Relatively  large 
values  were  recorded  at  the  850-mb  and  700-mb  levels  in  this  same  area. 
Again  the  rapidity  of  movement  of  the  front  and  associated  values  of  ad- 
vection and  short  length  of  the  data  period  prevented  detailed  investiga- 
tion but  suggests  the  possibility  of  a  meteorologically  significant  rela- 
tionship. 
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5.     Conclusions 

1.  Where  there  is  curvature  in  the  pressure  patterns  involved,  the  non- 
divergent  velocity  is  a  better  approximation  to  the  wind  than  the  geo- 
strophic  value, 

2.  No  significant  change  is  made  by  combining  the  non-divergent  and 
divergent  winds  for  the  uses  investigated. 

3.  A  relationship  possibly  exists  between  the  divergent  wind  and  front- 
al areas  that  could  be  utilized  in  locating  such  phenomena      This  con- 
clusion is  based  on  very  limited  data  but  offers  an  interesting  avenue 
for  further  investigation „ 
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FIGURE  1 
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FIGURE  2 
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FIGURE  12 
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APPENDIX  I 


Pillow  is  defined  as 


P=      Z(ft-B) 


n  (i-D 


X 

The  root-mean- square  differnece  (RMSD)  is  determined  as  follows 


HMSD- 


~\ 


y  =  n  =  3«)fc8 


In  equations  (1-1)  and  (1-2),  An  and  Bn  represent  the  point  values  of 
the  fields  being  compared.     In  this  example  the  A  value  is  a  geostrophic 
value  of  speed  or  direction  and  B  is  the  corresponding  value  of  either  non- 
divergent  or  total  wind. 
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